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Abstract 

We derive a semiclassical time evolution kernel and a trace formula for the Dirac 
equation. The classical trajectories that enter the expressions are determined by the 
dynamics of relativistic point particles. We carefully investigate the transport of the 
spin degrees of freedom along the trajectories which can be understood geometrically 
as parallel transport in a vector bundle with SU(2) holonomy. Furthermore, we give 
an interpretation in terms of a classical spin vector that is transported along the 
trajectories and whose dynamics, dictated by the equation of Thomas precession, 
gives rise to dynamical and geometric phases every orbit is weighted by. We also 
present an analogous approach to the Pauli equation which we analyse in two different 
limits. 
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1 Introduction 



Most semiclassical approaches to the Dirac equation so far aimed at an extension of the 
WKB method, with the expectation that some kind of Bohr-Sommerfeld quantisation con- 
ditions would emerge in relativistic quantum mechanics. The earliest such approach is 
due to Pauli Q, who succeeded in showing that the phase of a WKB spinor is given by 
a solution of the Hamilton- Jacobi equation for relativistic point particles. But he could 
determine the amplitude of the semiclassical spinor only in some special cases. Although 
the programme mentioned above has been very successful in nonrelativistic quantum me- 
chanics, where it leads to the so-called Einstein-Brillouin-Keller (EBK) quantisation [Q, 
establishing semiclassical quantisation conditions in the case of the Dirac equation was 
found to be obstructed by the occurence of certain phases. Investigating multicomponent 
wave equations Yabana and Horiuchi [H] noticed that geometric phases play an important 
role in this context and have to be incorporated in appropriate quantisation conditions. 
Using path integrals instead of the WKB method, Kuratsuji and Iida [|], ||] realised that 
an inclusion of the geometric phases in the symplectic form on classical phase space offers 
a possibility to arrive at quantisation conditions. A general theory of semiclassical quan- 
tisation for multicomponent wave equations in arbitrary dimensions, which derive from 
classical Hamiltonian matrices with no globally degenerate eigenvalues, was developed by 
Littlejohn and Flynn U [/]. This method, however, does not apply to the Dirac Hamil- 
tonian since the eigenvalues of the associated classical Hamiltonian matrix are twofold 
degenerate. Indeed, Emmrich and Weinstein [BJ found that in the degenerate case integra- 
bility of the corresponding classical dynamics is not a sufficiently strong condition to allow 
for an extension of EBK quantisation to the case of multicomponent wave equations. 

Moreover, even if successful, the procedure described above would not apply to systems 
whose classical limit is nonintegrable. In this paper we will therefore follow an alternative 
approach in that we investigate the semiclassical time evolution and then derive a semi- 
classical trace formula. This method is in the spirit of Gutzwiller's semiclassical treatment 
of the Schrodinger equation || (see also [[TTJ), in which the quantum mechanical density 
of states is set into relation to a sum over the periodic orbits of the corresponding classical 
system. The geometric phases mentioned above also appear in our approach in that they 
represent the spin transport along classical orbits in the trace formula. The advantage of a 
trace formula approach is not only that it is applicable to both integrable and chaotic sys- 
tems, but that it also provides the basis for efficient semiclassical quantisation conditions. 
In the context of quantum chaos extensive studies in this direction have been undertaken, 
see e. g. 0, Furthermore, semiclassical trace formulae are the primary tools for a 
semiclassical theory of spectral statistics, see e. g. [14], [15|, [16 . 



Before we go into more detail we want to take the opportunity to clarify our point 
of view regarding the semiclassical limit on which we base the following investigations. 
In general we will consider the mathematical limit h — > 0, with the understanding that 
in a given physical situation an equivalent limit in terms of physical quantities such as 
controlable external parameters has to be taken. The effect of such a limit can then be 
expressed in terms of the spectrum of the Dirac Hamiltonian as follows. The l.h.s. of the 
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trace formula that we are going to derive reads 



E n — E 

h 



where the E^s are the quantum mechanical eigenvalues, which depend on h, and E denotes 
a variable parameter. The smooth test function g decreases faster than any power for large 
arguments so that the main contribution to the sum comes from eigenvalues within an 
interval of length proportinal to h, 

E + hw x < E n < E + ftw 2 . (1.2) 

In the formal limit h — ■> the Weyl law forces the spectral density to increase in such a 
way that, although the length of the interval ( |1.2|) shrinks to zero, infinitely many terms 
contribute to In other words, the semiclassical limit corresponds to the limit of an 

increasing spectral density. The latter can possibly be achieved in a variety of different 
ways, which sometimes makes it necessary to vary several external parameters simultane- 
ously. In case the limit h — > is accompanied by further limits, one then has to ensure 
that asymptotic expansions are uniform in the quantities that are involved in the further 
limits. Throughout this paper we will understand the semiclassical limit as being involved 
with the formal asymptotics as h —>■ 0. In the last section, where we investigate the non- 
relativistic approximation of the results obtained for the Dirac equation, we briefly discuss 
an example with a second, simultaneous limit. 

In the following sections we will develop the steps that are necessary to derive a trace 
formula for the Dirac equation in some detail. Basically we follow the method introduced by 
Gutzwiller |T7], ITS), [5] in the case of the Schrodinger equation. It turns out that regarding the 
translational degrees of freedom Gutzwiller's approach can indeed be taken over. The novel 
features that arise in the case of the Dirac equation derive from the spin degrees of freedom 
and their coupling to the translational dynamics. In section |2] we first fix our notation 
and recall the basic properties of the Dirac equation relevant for the following. Then we 
review the general role of periodic orbits in semiclassical trace formulae and introduce a 
regularisation procedure (cf. [[19], which allows to obtain convergent trace formulae. We 
moreover recall how to cut off the essential spectrum of the Dirac Hamiltonian, which is 
present in most physically relevant situations and typically covers (— oo, —mc 2 ] U [mc 2 , oo). 

Section H is devoted to the derivation of a semiclassical time evolution kernel in the spirit 
of the Van Vleck formula known for the respective kernel for the Schrodinger equation. In 
the present context we find it convenient to represent the time evolution kernel in terms of 
an oscillatory integral. This method was developed for the study of scalar wave equations 
in the context of microlocal analysis, see e. g. [plfl, and subsequently found application to 



the development of several trace formulae [22, 23, [L9l E0|. In the case of the Schrodinger 



equation it leads to the same result as Gutzwiller's original derivation |17| which employed 
a stationary phase approximation of a Feynman path integral. The approach that we are 
going to follow is similar to the usual WKB method and results in equations that determine 
the coefficients of an /i-expansion of the time evolution kernel. The presence of spin is 
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reflected in these equations through their matrix character. To lowest order one obtains 
as a condition for their solvability two Hamilton- Jacobi equations, which correspond to 
the classical dynamics of relativistic point particles with positive and negative kinetic 
energy, respectively. The condition that arises in next-to-leading order in h is usually 
called transport equation. The latter can be reduced to two differential equations for 2 x 2 
matrices describing the transport of the spin degrees of freedom along particle orbits. The 
solutions of these spin transport equations as well as the solutions of the Hamilton- Jacobi 
equations finally determine the leading order of the semiclassical expansion for the time 
evolution operator. It is remarkable that the classical Hamiltonians do not include any term 
corresponding to forces acting on the magnetic moment of the spin. This fact was already 
realised by Pauli in his WKB treatment of the Dirac equation, see also a related discussion 
24J . For this reason Pauli's method was subsequently criticized by de Broglie p5[ , 



m 



who argued that electromagnetic moments linked with the spin were classical quantities 
and therefore should be present in a semiclassical approximation. This objection was later 
clarified by Rubinow and Keller . They pointed out that the moments of an electron are 



proportional to h so that in leading order as h — > the influence of spin on the trajectories 
vanishes. Furthermore, they showed how to obtain the equation describing the Thomas 
precession [27] of a classical spin, which is also known as the BMT-equation [28|, from the 
transport equation. Since it only contains the ratio of the magnetic moment and the spin, 
H cancels from the equation of Thomas precession which therefore provides the correct 
description of a classical spin. 

In section f| we analyse the spin transport in more detail, where we mainly focus on two 
aspects. We first discuss the geometric terms that accompany semiclassical asymptotics of 
multicomponent wave equations as they follow from the transport equation. Their structure 
is of a similar form as discovered by Littlejohn and Flynn J7| in the case of wave equations 
with classical Hamiltonian matrices that have no (globally) degenerate eigenvalues. In 
particular, one contribution is identified as being of the same type as the Berry phase 
|29| , |30| , |3T| appearing in adiabatic approximations. However, since in the case of the Dirac 
equation the eigenvalues of the classical Hamiltonian matrix are twofold degenerate, the 
U(l)-holonomy factors of Littlejohn and Flynn are replaced by corresponding SU(2) terms. 
Our results are found to be in accordance with the general discussion of the transport 
equation for multicomponent wave equations by Emmrich and Weinstein ||. In addition, 
these authors revealed the global geometric meaning of all terms that contribute to the 
total holonomy following from the transport equation. As a second point we discuss how 
to express also the spin contributions to the semiclassical time evolution kernel in terms 
of classical quantities. To this end we introduce a classical spin vector s as an expectation 
value of a time dependent spin operator. It then follows from the spin transport equation 
that s has to fulfill the classical equation of Thomas precession. Thus, up to a phase 
factor the desired solution of the spin transport equation is determined by the classical 
spin precession along a given particle orbit. We also show that the additional phase factor 
is composed of a dynamical part associated with the energy of a classical magnetic moment 
in given electromagnetic fields, and a geometric part which is of the type discovered by 
Aharonov and Anandan [R2 . 
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In section [| we derive the semiclassical trace formula for the Dirac equation, which is 
our central result, by Fourier transforming from the time domain to the energy domain 
and by subsequently taking the trace over spatial variables and matrix components. The 
periodic orbits that enter the trace formula are determined by the relativistic dynamics 
of classical point particles without internal degrees of freedom. The influence of the spin 
appears through two phase factors every periodic orbit is weighted by. One of these phases 
measures the change in the direction of the classical spin after this has been transported 
along a periodic orbit, whereas the second one contains the phase described above as being 
composed of a dynamical and a geometric part. 

Finally, in section ^ we consider the nonrelativistic limit of the results obtained in the 
previous sections. It turns out that the leading order as c — > oo coincides with the result 
of an application of the above programme to the Pauli equation. This equation is also 
of independent interest because in applications it is often used to investigate spin-orbit 
coupling. Based on the method developed in [7|], this important effect is e. g. treated semi- 
classically in |33| , |34| . Frisk and Guhr |34j] introduced a trace formula for a nonrelativistic 



Hamiltonian that describes spin-orbit coupling by modifying the Berry- Tabor trace formula 
p5| for classically integrable systems without spin. To be able to use the formalism of |7| 
they kept Bohr's magneton \i = fixed as h — ► 0; otherwise the corresponding classical 
Hamiltonian matrix would have a twofold degenerate eigenvalue. This procedure is in con- 
trast to our method which allows to treat H on the same footing in all terms. Fixing Bohr's 
magneton can be regarded as simultanously performing the limit h and the limit of 
an infinite coupling of spin to the translational degrees of freedom. In the semiclassical 
expressions the spin precession then decouples adiabatically from the translational motion. 
For the Pauli equation that describes a coupling of spin to an external magnetic field we 
compare two ways of performing the semiclassical limit: (i) h — > and (ii) h — > 0, \B\ — > 
with h\B\ = const. In the first case the nonrelativisitic limit of the result for the Dirac 
equation emerges, where the geometric terms in the spin transport are SU(2)-holonomy 
factors. In the second case the adiabatic decoupling of the spin motion results in geometric 
terms that yield a U(l) holonomy such that the Berry phase of a precessing spin is recov- 
ered. Both results being different implies that the semiclassical limit /I — ^ is not uniform 
in the magnetic field strength, a result that sheds some light on de Broglie's criticism of 
Pauli's approach. 



2 Semiclassical asymptotics and classical trajectories 

In this section we will present the basis for our subsequent discussion of semiclassical 
methods for the Dirac equation. In the following our focus will be on relativistic particles 
of charge e and mass m with spin 1/2 in external static electromagnetic fields. Thus the 
relevant Dirac equation reads 

ih—{x,t) = H D ^{x,t) , (2.1) 
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with the quantum Hamiltonian 



H D ■= ca ( - -A(x) ) + (3mc 2 + e tp(x) , (2.2) 



1 

which is a matrix-valued differential operator of first order. Here <p and A denote the 
electromagnetic potentials such that the corresponding fields are given by E(x) = — Vg<p(x) 
and B(x) = x A(x). The Dirac algebra is realised by the 4x4 matrices 

<j\ , „ ( 1 2 



S=| s Q j and | , (2.3) 



where a denotes the vector of Pauli matrices, and It. is the k x k unit matrix; see 36 for 



further details. The Hamiltonian (|2.2|) can be realised as 



H D = H D [-V,,x) , (2.4) 



i 



where 



Hd (p,x) := ca (j) A(x)j + f3mc 2 + e(p(x) (2.5) 



is the symbol matrix, in the sense of Weyl quantisation, of the operator H D (see e. g. |37j]). 
In case the potentials if and A satisfy suitable regularity conditions, the Dirac Hamiltonian 
H D , when defined on the domain C£°(IR 3 ) ® C 4 in the Hilbert space H := L 2 (1R 3 ) ® C 4 , 
is essentially self-adjoint; see e. g. |3£| ^9|. In the following we will always deal with its 
self-adjoint extension which we also denote by Hp. 

Solutions *$?(x, t) of the Dirac equation ( j2.1|) , with initial conditions *&(x, 0) = ^o(^) £ 
TC, are obtained as U(t)^o — e~^ Hot ^Q. The time evolution operator U(t) can then be 
represented by its Schwartz kernel K(x,y,t), so that 

t)= f K{x, y, t) * {y) d 3 y . (2.6) 

This matrix-valued kernel is obviously required to solve the Dirac equation, with initial 
condition 

Km K(x, y, t) = 1 4 5{x - y) . (2.7) 

Our first major goal, to be dealt with in the next section, will be to derive a semiclas- 
sical representation of the time evolution kernel K (x, y, t) in the spirit of the Van Vleck 
formula for the respective kernel of the Schrodinger equation. In the latter case such a 
representation is usually derived from a Feynman path integral, to which the method of 



stationary phase is applied [17|. Here we prefer an alternative semiclassical approach that 



makes use of a representation of the kernel in terms of an oscillatory integral. This method 
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has previously also proven useful in the case of the Schrodinger equation, see e. g. HO, EL6 



in which it served as a basis for a mathematically rigorous proof |T9|, of the Gutzwiller 
trace formula 



m. 



We therefore now introduce the matrix-valued oscillatory integral 
1 



K(x, y, t) 



real-vf 

4x4 matrices with semiclassical expansions 



a+(x, y, t; £) e^ + (^ + a^(x, y, t; e** _ ^« d 3 £ , (2.8) 



where ± are real-valued smooth phase functions that are independent of h, and are 



a%(x,y. 



oo , ^ \ k 

«) = E t a±(x,y,t;0+O(h c 



[2.9) 



According to the general theory of Schwartz kernels, K(x, y, t) is a distribution kernel so 
that the integral ( |2.8| ) has to be interpreted in a distributional sense, see e. g. [^DJ for 
details in the scalar case. In order to account for the initial condition (12.71) one chooses 



± (^^,O;O = {x-y)£, , 



and 



at{x,y,0-^) + a k (x,y,0;£) 



U if k = , 
if k > 1 . 



(2.10) 



(2.11) 



Since the kernel ( |2.8| ) has to solve the Dirac equation, one obtains conditions for the phases 
± and the coefficients appearing in the matrix- valued amplitudes a^. A detailed 
discussion of these equations will be postponed to the next section. Here we only remark 
that the phase functions have to satisfy the Hamilton- Jacobi equations 



H± (V^ ± (x, y, t; 0, i) + ^(x, y, t; () = 
with the classical Hamiltonians 



H (p, x) = etp(x) ± \ c 2 [p A(x) ) + m 2 c 



(2.12) 



(2.13) 



of relativistic particles with positive (+) and negative (— ) kinetic energies, respectively. As 
can readily be verified, the functions H ± (p, x) are the two, twofold degenerate, eigenvalues 
of the symbol matrix ( |2.5| ). A posteriori, the occurrence of two Hamilton- Jacobi equations 
( p,12| ) justifies the choice (|2.8|) of the oscillatory integral with two additive contributions. 

Due to the form fl2.12|) of the Hamilton- Jacobi equations and the initial conditions 
( p.lO| ), the variable y can be separated according to 



<P ± (x,y,t;0=S ± (x,£,t) 



(2.14) 
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so that the functions S ± (x,!;,t) also solve the Hamilton- Jacobi equations ( |2.12| ), but with 
initial conditions S' ± (x, £, 0) = x^. From general Hamilton- Jacobi theory, see e. g. |4l], |0| , 



it is known that therefore S (x, £, t) are generating functions for canonical transformations 
(p, x) i— > (£,-?). For ease of notation we restrict the following discussion to the index +. 
Then (p, x) are the end points of the solution (P(t'),X(t')) of Hamilton's equations of 
motion, generated by the Hamiltonian H + , with initial condition (£, z). This means that 

(P(0),X(0)) = (£>) and (P(t),X(t)) = (j),x) . (2.15) 
In order to explicitly incorporate the initial conditions we will also use the notation 

P(t') = P{t'- £ z) and X(t') = X(t'; £ z) for < t' < t . (2.16) 

The fact that S + (x, £, t) is a generating function for the associated canonical transformation 
moreover implies the relations 

p = (x, £ t) and z = V^S + (x, £ t) . (2.17) 



For sufficiently small times, \t\ < t C} the Hamilton- Jacobi equations ( |2.12j ) are known to 



possess unique solutions S ± (x, £, t), see e. g. ipLOf . However, it is well known that at some 



critical time t c a caustic may arise so that for \t\ > t c the solutions of the Hamilton- 
Jacobi equations are no longer unique. The representation ( [2.8|) of the time evolution 
kernel for \t\ > t c can then only be a local one. A global object has to be constructed 
by gluing appropriate local representations together. This procedure requires consistency 
conditions, which are reflected by the presence of Morse indices in the final semiclassical 
expression to be derived in section [| For the Schrodinger equation this was already noticed 



by Gutzwiller |[L7| . A mathematically rigorous construction employing the Maslov bundle 
can be found in fL9f . For simplicity, we will neglect this problem below in that we continue 
to work with (|2.8|) and only in the end introduce the appropriate phase factors. 

Since the principal applications of semiclassical trace formulae deal with eigenvalues 
of quantum Hamiltonians, in particular with their semiclassical determination and their 
statistical properties, respectively, one would like to isolate the point spectrum of Hp from 
its essential spectrum. In case the potentials <p and A vanish towards spatial infinity, 
or under suitably weakened conditions, the essential spectrum is known to be given by 
IR\(— mc 2 , mc 2 ), see e. g. |38|, so that for the following we assume Hn to have, possibly 



infinitely many, eigenvalues of finite multiplicities in the interval (—mc 2 ,mc 2 ). Given 
any interval / = (E a , Ef,) that contains only isolated eigenvalues, we introduce a smooth 
function x(E) that is supported in J, i. e., it vanishes outside of J, such that x{E) = 1 
on a suitably large subinterval. In particular, if there is no accumulation of eigenvalues 
at E a or Et>, one can achieve that x(-^n) — 1 for all eigenvalues E n of Hp. Otherwise, as 
e. g. for the Dirac Hamiltonian of the hydrogen atom, one can enlarge the support of \ 
arbitrarily towards ±mc 2 so that arbitrarily many eigenvalues can be taken into account. 
Then the operator x{Hd), defined by the functional calculus given by the spectral theorem, 
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is bounded and self-adjoint; see e. g. [4(J for the scalar case. This operator has a purely 
discrete spectrum with eigenvalues x(E n ). Therefore the Schwartz kernel K x (x, y, t) of the 
truncated time evolution operator U x (t) := x{Hd) U(t) has a spectral representation 



K x (x,y,t) = J2x(E n )MxM(y)e~* Ent 



[2.18) 



in terms of the orthonormal eigenspinors \l/ n of Hp, where ^ denotes the hermitian adjoint 
of \l/ n . According to the definition of the truncated time evolution operator, which implies 
that 



the truncated kernel is obtained from the non-truncated one as 

K x {x, y, t) = x(H D ) K(x, y, t) 
1 



X(H D ) a+(x, y, t; £) e^ + ^« + a^(x, y, f ; e**"<^ 



(2.19) 
(2.20) 

d 3 £; 



here x{Hd) always acts on functions of x. In the framework of Weyl calculus one can iden- 
tify the symbol of the operator x{Hd) to possess a semi classical expansion with principal 
symbol (leading term) x(H D (p,x)), see |40|| . Thus, an explicit calculation yields that 



M ( x (H D )ae^) (f) = x (h d {V s <j>{x), x)) a{x)+0{h) 



(2.21) 

for any sufficiently regular matrix-valued function a and real-valued function <p. This 
calculation is closely parallel to the respective result in the scalar case, which can be found 
in 



2l|] . Notice that in |2l] a different quantisation is used which, however, yields to lowest 



semiclassical order the same result as Weyl quantisation. Thus, to leading order in h, the 
truncated time evolution kernel reads 



K x (x,y,t) 



X (H D (V^ + ,x))a+ (l+0(h))ei+ 
+ X {H D {Vz<j>-, x)) a (1 + 0(h)) e^' 



(2.22) 



d 3 e 



In order to prepare for the semiclassical trace formula to be dealt with in section |5| 
we now introduce a regularisation of the truncated time evolution operator. To this end 
consider a smooth test function q G C°°(R) such that its Fourier transform 



g(t) 



g(E) e^* dE 



(2.23) 



is smooth and compactly supported. Then in particular, both g and g are test functions 
from the Schwartz space iS(M). We now define the (bounded) operator 



U x [g) :-- 



1 

2^ 



+oo 



g(t) U x (t) dt 



(2.24) 
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whose trace can be calculated with the spectral representation ( |2.18| ) of the truncated time 
evolution kernel as 



Tr UM = ^J R J e(t) trK x (x,x,t) dtd 3 x = J^xiEn) g 



(2.25) 



Here Tr means the operator trace on L 2 (R 3 ) ® C 4 , which includes a trace over the matrix 
components; the latter is denoted by tr. The linear map g i— > TrJ7 x [^] defines a tempered 
distribution, denoted by Trf/ X [-] G «S'(M), if the sum on the r.h.s. of ( j2.25| ) converges 
absolutely. To this end one requires that 



x(E n ) g ( f 



O 



n 



l+e 



e > 



n 



oo . 



(2.26) 



The simplest case of finitely many eigenvalues in the support of x obviously poses no 
problem. If, however, eigenvalues accumulate at some E acc £ [—mc 2 ,mc 2 ], the truncation 
X has to be chosen such that x{E aC c) — and, moreover, x{E) vanishes sufficiently fast as 
E — > E acc in order to fulfill ( |2.26| ). We now evaluate the distribution Tr[/" x [-] on the test 
function g(t) e~% Et and thus obtain the relation 



E n — E 

h 



l 

2^ 



g(t)e^ Et trK x (x,x,t) dtd^x 



(2.27) 



The semiclassical trace formula we are aiming at results, if for the truncated time evolution 
kernel in ( [2.27D one introduces a semiclassical representation and calculates the integrals 
with the method of stationary phase. The details of this procedure will be carried out in 
section [5]. Here we only remark that upon introducing the representation ( [2.22 ) for the 
kernel, in leading semiclassical order one has to compute the integrals 



g(t) X (H D ($3<f& x, t; |5, £)) 4{x , x, t; f) e ^At^+Et] df d 3 x ^ 

(2.28) 



According to the method of stationary phase, see e. g. [ED, ET] , all contributions to (|2.28| 



that exceed 0{ff°) as h — > are determined by the stationary points, in the variables 
(£, x, t), of the phase ± (x, x, t;£) + Et. These stationary points are solutions of the 
equations 



V^{x , y, t; + V^x, y, t; 
(x , x , t; I) + E = 







J y=x 



dt 



(2.29) 



Vxf ) ± (x,x } t;O = . 
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If one now recalls the connection ( [2.141 ) of the phase ± to the generating function S ± , one 
obtains the equivalent conditions 

V^ ± (f,e,t)=x , VtS ± (x,U)=Z , -^-(£,Z,t)=-E , (2.30) 

to be fulfilled by stationary points (£, x, t). A comparison with the relations ( |2.17| ) now 
yields the conditions £ = p and x = z, so that the stationary points determine periodic 
solutions, with energy E, of Hamilton's equations of motion generated by the Hamiltonians 
( |2.13| ). We denote these periodic orbits by 7?. Their periods, given by the t-components 
of the corresponding stationary points, are called T 7 ± . Notice that our requirement on the 
test function g to have compact support implies that only periodic orbits up to some finite 
period, \T^± \ < T max , contribute. However, since the support of g can be made arbitrarily 
large one can manage to include as many periodic orbits as desired. As our first observation 
on the way towards a semi classical trace formula we thus now conclude that, apart from 
terms of 0(h°°), in the semiclassical limit h — > all contributions to the l.h.s. of ( 2.27] ) are 
due to classical periodic orbits of energy E of relativistic point particles. 



A further observation can be made with (|2.28|) if the set of stationary points (£, x, t) 
divides into a sequence k G Mo, of m^-dimensional smooth submanifolds of IR 3 x IR 3 x R 
such that the phase S (x, £, £) — x£ + Et is non-degenerate transversal to the manifolds 
Mfc. The latter condition means that the matrix of second derivatives of the phase with 
respect to (£, x, t) has constant rank 7 — irtk when restricted to Mk- The flow generated by 
if ± is then called clean, see |L9], for further details. In this situation the set of periods 



is known to be discrete in IR [ 52 ] . Thus the periods cannot accumulate at some finite value 
T. Under these conditions the method of stationary phase implies that each manifold Mk 
yields a contribution 

00 

h^AM^e^+Oih 00 ) , with A Mh {h) = ^2n?A MhJ , (2.31) 



to (|2.27D . Here Su k is the action of any periodic orbit 7^ contained in Mk, which can 



be computed as the integral of the canonical one- form pdx along any closed path on Mk- 
The coefficients AM k ,j are independent of ft, and arise in the method of stationary phase as 
certain integrals over Mk- In the case of an isolated periodic orbit 7^ the manifold Mk is 
given by the points on the primitive periodic orbit associated with 7^, and thus rrik = 1. 
If the phase is non-degenerate transversal to an isolated periodic orbit, the latter is called 
non-degenerate. For this situation we will explicitly calculate AM k ,o hi section [5[ Another 
case that can be dealt with explicitly concerns the hypersurfaces 

fi± :={(££); H ± (p,x) = E] (2.32) 

of constant energy E in phase space. The points (p, x) G are obviously periodic 
under the flows generated by H , with trivial period To = 0. Since according to the 
above non-degeneracy condition we assume that M$ = x {0} are (five dimensional) 
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smooth manifolds, E is required to be a regular value for the respective flows. Moreover, 
the associated leading terms in ( |2.31|) are of the order h~ 2 . The explicit calculation of 
A M ± o will also be performed in section |5|. As a final remark let us mention that in case 
the dynamics generated by H ± are integrable, phase space foliates into three dimensional 
invariant tori such that the respective manifolds My. are given by the rational ones among 
these tori. Thus = 3 so that the leading terms in (|2.31|) are of the order h~ l . 



3 Semiclassical time evolution 



In this section we will determine semiclassical representations for the time evolution kernel 
K(x,y,t) and its truncated version K x (x,y,t), respectively. In a first step we derive 
necessary conditions that must be imposed on the matrix-valued amplitudes and on 



the real-valued phases ^ in order that the oscillatory integral ( |2.8|) be a possible ansatz 
for the time evolution kernel. In a second step we then employ the method of stationary 
phase to (|2.8|) , and from the result we obtain a semiclassical expression in the spirit of the 
Van Vleck formula for the respective kernel of the Schrodinger equation. 

The requirement that (|2.8|) , together with the expansion ( |2.9|) , be a semiclassical repre- 
sentation of the kernel K(x, y, t) to all orders in h means that the oscillatory integral shall 
fulfill the Dirac equation (|2.1|) up to terms 0(h°°). We therefore act with Hjj — ihjjr on (|2 



after having introduced the expansion (|2.9| ), and then group terms of like orders in h. The 
phases ^ and coefficients af now have to satisfy equations that result from demanding 
that all coefficients of powers h k , k = 0, 1, 2, . . . , vanish. To lowest order {k = 0) one thus 
obtains the equation 



Hjj{ Yj:< > + ( x. ff. /-:<).,?)+ f^L (f ? £ f, f) 



H D ( V £ (x, y,t;£),x) + -^-(^ V, 



(3.1; 



in which Hd(p,x) denotes the symbol matrix (|2.5| ) of the quantum Hamiltonian Hp. For 
the following it turns out to be convenient to demand an individual vanishing of the terms 
with index + and — , respectively. Since the two twofold degenerate eigenvalues of the 
hermitian 4x4 symbol matrix Hd(p,x) are given by H ± (p,x), see ( |2.13| ), the condition 
( j37TD can be fulfilled as soon as the matrices Oq are suitably composed of eigenvectors of 
Hd^x^jX). Upon diagonalising the symbol matrix Hd(p,x) one obtains an orthonor- 
mal basis for C 4 (endowed with the canonical scalar product) that consists of eigenvec- 
tors {ei(p, x), e 2 (p, x)} with eigenvalue H + (p,x), and eigenvectors {fi(p,x),f 2 (p,x)} with 
eigenvalue H~(p,x). We now define the 4x2 matrices V(p, x) and W(p,x) whose two 
coloumns are given by the vectors e\ (p, x) , e 2 (p, x) and fi(p,x),f 2 (p,x), respectively. In 
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explicit terms these matrices read 

I f[e(p,x) +mc 2 ]tr 

V{p,x) 
W(p, x) 

where 



y/2e{p,x)[e{p,x) +mc 2 } \ W~ eA{x)\ a 
1 ( [cp — eA(x)]a 

y/2e(p,x)[e(p,x) +mc 2 } \~Hp, x) + mc 2 ] 1 2 



(3.2) 



e(p,x):=^(cp-eA(x)) + m 2 c 4 = H + (p, x) — e<p(x) . (3.3) 

Since the eigenvectors {ei, e 2 , fi, / 2 } are chosen to be orthonormal and form a basis for C 4 , 
one obtains the relations 

\rtV = ± 2 , W ] W = l 2l VV ] + WW ] = t A . (3.4) 

Moreover, the projectors P±(p,x) onto the eigenspaces corresponding to the eigenvalues 
H (p,x), respectively, are given by P + = VV* and P_ = WW^. 

Following Q2.14j ), we now choose the generating functions S^^ to replace the phase 
functions <p ± . If one then introduces suitable 2x4 matrices V and W, 

4(x, y, t; i) = V (v £ S + (x, £ t), f) V(x , £ 1) , 

a ~ (x , y, t; £) = W (v g S~ (x, £ 1) , f) W(f , £ t) , (3.5) 

the condition ( |3.1| ) is fulfilled as soon as S' ± , and hence also the phases ± , satisfy the 
Hamilton- Jacobi equations 

F± (V^Or,^),^ + -gf&^t) = • (3-6) 

In the definition of the matrices V and we anticipated the fact that the coefficients 
are independent of y, see the discussion below ( |3.8|) . At the moment an explicit expression 



for V and W is not needed. After having applied the method of stationary phase to (|2.8| ) we 
will specify them further. Here we only remark that in order to fulfill the initial condition 
( [2. 11| ) we demand that 

V(x, £ 0) = 0(£ x) and W(x, £ 0) = W\£, x) , (3.7) 

so that <2q (x, y, 0; £) = -P± (£,:?). In conclusion, the condition ( |3.1| ) to lowest order in h 
requires the phases (f)^ to solve the Hamilton- Jacobi equations ( p. 12 ). Due to the initial 



condition ( |2.10| ) these functions are therefore now completely fixed. 

Our next goal is to determine the leading contributions to the amplitudes. The 
condition ( p.l| ) appearing to lowest order in h only requires these coefficients to be of a 
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certain structure, see ( |3.5| ). They are, however, completely fixed by the conditions imposed 
in next-to-leading order. Regarding all higher orders h k , k = 1, 2, 3, . . . , one observes that 
the equations obtained from the procedure described before Q3.1|) can be given in a uniform 
manner. They read 



H D (V s S + ,x) + 
H D (V s S~,x) + 



dS+ 
~df 
dS_ 
dt 



9 

ca\7 s + — 
at 

d 

ca\7 s + — 
at 



{3.1 



Since these equations, as well as the initial conditions ( |2.11| ), are independent of y, the 
coefficients do not depend on y either. Beginning with k = 1, the hierarchy (|3.8| ) of 
equations can principally be solved order by order. To this end, for each k we multiply 
the two equations with V A ^(V^S ,+ ,x) and W^(WgS~ , x), respectively, from the left. Since 
V^H D = V^H + and W^H D = W^H', the Hamilton- Jacobi equations Q3.6p then imply 
that 



V\V s S + ,x) 
W\V £ S~,x) 



c + — 



c aVi; + — 



d_ 
at 

o_ 

dt 



Ofc_ 1 (x,y,t;|) = . 



(3.9) 



If one started with as given in (|3.5| ), one now could in principle determine all higher 
coefficients from (|3.8| ) and ( |3.9|) recursively. Our ambition is, however, limited to obtain 
a semiclassical expression for the time evolution kernel to leading order. In the following 
we therefore restrict our attention to the case k — 1 and, moreover, only present the case 
with index + explicitly. 

Expressing as indicated in (|3.5| ), the l.h.s. of the equation (|3.9|) for k — 1 and index 
+ can be viewed as an application of the matrix-valued differential operator 



x) 



d_ 

dt 



V(V s S + ,x) 



(3.10) 



to the matrix- valued function V(x, £, t). A direct calculation shows that ( p.lOj ) can be 
expressed as 



VpH+(V s S + , x)V 3 + j t + [v f H + (V s S + , f)] + iM+(V g S + , x) 
where we defined the traceless hermitian 2x2 matrix M + as 



(3.11) 



M + (p,x) :-- 



ec 



2e(p,x) 



B(x) + ° (e(x) x(p- -A(x) 

e(p, x) + mc 2 \ \ c 



a . 



(3.12) 
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We remark that since M + is a linear combination of the Pauli matrices, \M + is an element 
of the Lie algebra su(2). In order to emphasize this point, and for later purposes, we 
introduce the following notation, 

1 



M+(p,x) 



R(p, x) a . 



(3.13) 



According to the usual convention, the resulting equation for V is called (first order) 
transport equation. At this point a comparison with the case of the Schrodinger equation 



seems to be instructive, see e. g. POJ. There the lowest order amplitude Qq is a scalar 



function that is required to be a solution of a transport equation which arises upon acting 



on a,Q Ch ^ with the equivalent to (|3.11[) , but where M + = 0. The solution, fixed by the 



initial condition ciQ Chr ' 



1 at t = 0, is well known to be 



,Schr 



{x, y, t; 



(3.14) 



Returning to the present situation described by ( ft.llp , the following ansatz for the solution 
of the transport equation therefore seems to be appropriate, 



V(x,U) 



'det 



/ d 2 S+ 
\dx k d£i 



(x,£t) U + {x,lt) . 



(3.15) 



As a consequence, the transport equation for V implies that the matrix-valued function 
U + has to solve the equation 



V^+(V^+, x)V s + ^ + iM+{V 3 3 + , x) 



(3.16) 



with initial condition U + (x, 0) = V^(^,x). 

It is not necessary to solve ( ^.16| ) in full generality because the amplitude a~Q enters 
the semiclassical expression for the time evolution only at stationary points £ 3 - of the phase 
<p + . We thus first employ the method of stationary phase to the oscillatory integral ( |2.8| ), 
in which we only take the lowest order contributions to the amplitudes into account. 
According to ( p.l4f j- fl2.17| ), the stationary points of the phase, i. e., the solutions £j of 



V^ix , y, t- Q = VgS^x, & t) - y = , 



(3.17) 



determine solutions {P{t'\ £j, y),X(t'; y)), < t' < t, of the classical equations of motion 
from y) to (p, x) in time t. Their projections X(t'; £j, y) to configuration space will also 
be denoted as 7^. At a stationary point the phase is given by Hamilton's principal function 
for the trajectory 7^" corresponding to £j, 



S ± (x,^,t) - y£j = J^L± (X{t')j{t')) dt' =: R 



"± (x, y, t) 

Ixy 



(3.18) 
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where L is the Lagrangian associated with H . 

The method of stationary phase requires the amplitudes ag , and hence U±, evaluated 
at the stationary points £j. In this situation x has to be considered as the end point of the 



corresponding trajectory 7. 



xy. 



1. e., x 



X(t). In ( p,16| ) the first two terms can therefore be 
understood as a total time derivative along X(t), 

%H + (P(t),X(t))Vg + - = X{t)V,+ - =: - . (3.19) 

In order to solve ( P-16|) at a stationary point £j we now introduce the ansatz 

[/ + (X(t),g,t) = d + (P(t),X(t))Vt(g,$ , (3.20) 

which immediately implies the initial condition 

d+(P(0),X(0))=l 2 (3.21) 

for the 2x2 matrix d + . Restricting ( |3.16| ) to the trajectory X(t) and using the abbreviation 
( P.19| ) then shows that this matrix is required to solve 



dt 



+ iM + (P(t),X(t)) 



d+(P(t),X{t)) = 



(3.22) 



along the trajectory. In the next section we will demonstrate that d + can be interpreted 
as a semiclassical propagator for the spin degrees of freedom. Occasionally, we will thus 
refer to ( ft.22|) as the spin transport equation. A formal solution of this equation can be 
given in terms of a time-ordered exponential, 



d+(t) = T exp <^ -i / M+(t') dt' 



(3.23) 



Since according to ( [3.12D \M + takes values in the Lie algebra su(2), the solution d + is an 
element of the group SU(2). Together with its connection to a classical spin, a further 
geometric interpretation of d + will be provided in the next section. 

Combining (p^5|), ( p,15| ), and ( |3.20| ) finally yields the lowest-order amplitude at a sta- 
tionary point as 



'det 



(■ 



d 2 S- 



(3.24) 



\dx k d£i 

V ( (x , g , t) , x) d + (V £ S + (x , g, t ) , x) (g , y) . 

for which the 2x2 matrix M_ 



This expression has an obvious analogue for the index 
is given by 



M_(p,f) :-- 



ec 



2e(p,x) 



B(x) - C - (E(x) x (p- e -A{x) 

e{p, x) + mc z V V c 



a 



(3.25) 
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Comparing ( |3.24| ) with ( |3.14| ) reveals that, apart from the occurrence of two classical Hamil- 
tonians, the only difference to the case of the semiclassical propagator for the Schrodinger 
equation is given by the appearence of the last three factors on the r.h.s. of ( |3.24| ), which 
we also abbreviate as Vtd + Vn. With the further short-hand 



Ixy 



d 2 ^ 



det 



Ixy 



dx k dyi 



x, y, t) 



1/2 



(3.26) 



the result of the method of stationary phase applied to ( |2.8| ) finally reads (for t 7^ 0) 



K{x, y, t) 



Ixy 



fly ' ySv [i + o(h)] 
+ J2 Wt d - W o D ix V e jR ^y~ l ^ y [1 + 0(h)] 



(3.27) 



Ixy 



where, as in the case of the Schrodinger equation, v+ denotes the Morse index of the 

; Ixy 

trajectory 7*. Notice that the two sums extend over the classical trajectories of relativistic 
point particles. The spin degrees of freedom, which manifest themselves in the Dirac 
equation through the matrix character of the Hamiltonian, enter in leading semiclassical 
order only through the factors V t d + V^ and W t gLWq , respectively. Since these contain no 
h, a classical interpretation of them seems to be in order. In the next section we will indeed 
obtain from ( p.22|) an evolution equation for a classical spin which is transported through 
the external fields along the particle trajectories 7^. 

As our final step in this section, we now modify ( |3.27j) to yield a semiclassical rep- 
resentation for the truncated time evolution kernel K x . According to fl2.22p , to this end 
one must apply the 4x4 matrix x{H D (W $S ± , x)) to the ampliude a^(x, y,t;£). Since the 
spectral representation Ho = H + P + + H~P^ of the symbol matrix implies 



X (H D ) = X (H + ) P + + X (H-) P = X {H + ) VV* + x(H-) W 
representing the amplitudes as in (|3.5|) immediately shows that 



x(H D )4 = x{H-- 



+ 



(3.28) 



(3.29) 



Since therefore the matrix-valued amplitudes are only multiplied by a scalar factor, which 
is moreover constant along any classical trajectory 7^, the only modification of (|3.27 ) 
consists of an inclusion of the factors x(-^ 7 ± ) under the sums over the trajectories, where 

E^± = H^Pty), X(t)) denotes the (constant) energy of 7, 



± 
xy 



17 



4 Geometry of semiclassical spin transport 



The semiclassical representation for the time evolution kernel that was derived in the previ- 
ous section is principally determined by the classical dynamics of relativistic point particles. 
This is true to the extent that it essentially suffices to solve the equations of motion gener- 
ated by the scalar Hamiltonians H^. Indeed, the semiclassically dominating /i-dependent 
phases occurring in ( |3.27| ) are completely fixed by Hamilton's principal functions i? ± of the 
translational motion. The latter decouples from the spin dynamics that, moreover, only 
contributes to the ^-independent amplitudes in (|3.27|) . 

In this section we will investigate the dynamics of the spin degrees of freedom more 
closely. In particluar, we will interpret the SU(2) matrix d + as a semiclassical time evolution 
operator for spin and, furthermore, relate it to a classical spin evolving according to a 
classical equation of motion. In order to prepare for this, we first have to clarify the 
geometrical setting of spin transport along classical particle trajectories further. To this 
end we consider the if + -eigenbundle E + over phase space. It is defined as the disjoint 
collection of the eigenspaces 



E 



(4.1) 



of the symbol matrix Hd(p,x) corresponding to the eigenvalue H + (p,x) for each point 
(p, x) in phase space. Since Hp is hermitian, the decomposition of C 4 into the H + - and 
if "-eigenspaces is orthogonal with respect to the canonical scalar product. This can hence 
be projected to the fibres so that E + is a C 2 vector bundle with structure group U(2). 

In terms of the orthonormal basis {ei(p, x) , e2{p, x)} used to define the matrix V(p, x) in 
(EO), a section in the eigenbundle can be represented as 



>(P> &) = ^2 v k(p, x) e k (p, x) = V(p, 



k=l 



(4.2) 



To lowest semiclassical order the spin dynamics is governed by the transport equation 
defined with the differential operator ( |3.10|) . This requires sections in the eigenbundle 
E + that are lifts of trajectories in phase space. We therefore now consider the following 
differentiation on sections v, which we evaluate at (p,x) = (V^S ,+ (x, £, t), x). Then 



v(VgS+,x) 



P+(V^ + ,f) 



d_ 

Of 



V(VzS + ,x)V\V s S + ,x) 



v(V s S+,x 
d 



(4.3) 



ca + 



dt 



V(V s S + ,x) 



vi(VgS + ,x) 
v 2 (V,S + ,x) 



is again a section in E + in the above sense. According to (|3.10| ) and ( p.ll| ) an expansion 
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of v in the gliding basis {ei, £2} yields the coefficients 



V 1 (V j ;S + ,X / 



+ iM+{V s S + ,x] 



v 2 (y s s + ,x) 



(4.4) 



Furthermore, a separation analogous to ( |3.15|) 



v(V s S + ,x) = Wdet 



/ d 2 S+ 
\dx k d£i 



(4.5) 



leads via (|3.16|) to the definition of a covariant differentiation T> + on the section u, 
V + u(V £ S + ,x) : = 



(4.6) 



If one now restricts the variable £ to the stationary points £j arising in (|3.17 ), the first 
two terms in ( |4.6|) again yield a differentiation along the trajectory 7+ y associated with £j, 
compare ( |3.22 ). One then also recognizes iM + as the su(2)-valued connection coefficient 
arising for the connection associated with T> + . A section u given along a solution of the 
classical equations of motion is therefore parallel, if it is a solution of 



V + u(P(t),X(t)) = with u k (P(0),X(0)) = u kfi . 
According to ( |3.22j ) such a solution can also be represented as 



Ul (P{t),X(t)) 
u 2 (P(t),X(t)) 



d + (P(t),X(t)) 



Ul,0 
«2,0 



(4.7) 



(4.8) 



Thus, geometrically d+ G SU(2) describes the parallel transport in E + defined by the 
connection arising from T> + . Along periodic orbits in phase space, d + therefore yields the 
holonomy of this connection. This geometric interpretation ensures that the combination 
V t d + Vq , which appears in the semiclassical expression ( j3.27|) for the time evolution kernel, 
is invariant under unitary base changes in E + . In physical terms, d + can be interpreted as 
the semiclassical time evolution operator for two-spinors in the representation defined via 
the gliding basis {ei, €2}- The unitarity of d + then implies that the norm |«i,o| 2 + |^2,o| 2 of 
the initial two-spinor (ui t o, ^2,o) T £ C 2 is preserved under this evolution. In the following 
we will therefore always consider normalised sections u in E + . 

We would now like to compare the above construction with the connections that appear 
in the analysis of Littlejohn and Flynn 0, and Emmrich and Weinstein 0. Since the 
eigenbundle E + is a (non-trivial) subbundle of the trivial C 4 bundle over phase space, a 
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natural connection in E + arises by projecting the trivial covariant differentiation of the C 4 
bundle to E + . Along a trajectory (P(t), X(t)) this construction reads 



P + (P(t),X(t)) 



dt 



u(P(t),X(t)) 



:V B u(P(t),X(t)) , 



(4.9) 



when applied to a section u in E + . In analogy to ( }4.6|) one then obtains a covariant 
differentiation in terms of the coefficients Uk with respect to the gliding basis {ex, 62}- The 
connection coefficient iM# that replaces iM + can be calculated as 



M B (p,x) 



ec 



2 e(p, x) [e(p, x) + mc 2 



a 



+ 



ec° 



2 e 2 (p, x) [s(p, x) + mc 



(p-^A(x))xE(x 

(p- lA{x)) x((p- lA{x)) x B{x)) 



(4.10) 



a 



That way one defines a connection on E + that bears some similarities to the adiabatic 
connection in quantum mechanics. The latter has been identified by Simon BO to produce 
the Berry phase [p9[] , see also ||. On this ground, Littlejohn and Flynn 0, 0] introduced the 
notion of a Berry term for the analogue to our Mb in the case of principal symbol matrices 
with no globally degenerate eigenvalues. In the present context, the twofold degenerate 
eigenvalue H + forces Mb to take values in su(2), which hence leads to a connection with 
SU(2) holonomy. We emphasize, however, that no adiabatic approximation is made after 
one has arrived at the transport equations ( |3.8|) and (|3.16|) , respectively. Nevertheless, 
following |6|, in spirit, we refer to Mb as the SU(2)-Berry term, although regarding 
adiabatic approximations this notation is sligthly misleading. 

The difference between the two covariant differentiations T> + and T> B can be expressed 
in terms of a connection coefficient Mc = M + — Mb- This expression, which Littlejohn 
and Flynn refer to as the no-name term, has been identified by Emmrich and Weinstein 
H to be related to the curvature associated with V B . They moreover showed that 



M c (p, x) = - l -V\p, x) {P + , H D - H + l 4 }(p, x) V(p, x) 



(4.11) 



where the Poisson bracket for two matrix- valued functions A, B on phase space is defined 



as 



{A, B}{p, x) := V p A(p, x) VxB(p, x) - V f B{p, x) V s A{p, x) . 



(4.12) 



Notice that the ordering of the matrices is important. Upon explicitly calculating the r.h.s. 
of ( P~np one can verify that M+ = M B + M c . 

We now define a time-dependent spin operator whose components act on C 2 , as 



E(t) :=d{(P(t),X(t))ad + (P(t),X(t)) , 



(4.13) 
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so that S(0) = a. According to ( |3.22| ), the dynamics of the spin operator E(i) is governed 
by 

±Z(t)=idl(t)[M + (t)cr-CTM + (t)}d + (t) . (4.14) 



Introducing R(p,x) as in ( |3.13| ), the evolution equation for E can be brought into the 
convenient form 

^E(t) = R(t) x E(t) , (4.15) 

which describes the precession of the vector E about the instantaneous axis defined by R. 
In order to obtain an object that can be considered as a classical spin, one introduces the 
expectation value of E in a two-spinor state prescribed by (wi,o, W2,o) T G C 2 , 



<f(i) := (^,0,^,0) X(t) M = (Mt)Mt))*M£) = I 2 Im(^(t) W2 (t)) I , (4.16) 



'2 Re(wi(t)u 2 (t))' 
Ul{t) u 2 (t)) 
\ Ul (t)\ 2 - \u 2 (t)\ 2 



which obviously also solves ( [4. 15[ ) , i. e. 



^ ?x (^-HFT^)( p "-^) xJ ) ' (417) 

where all quantities have to be taken along a given trajectory 7+ This classical equation 
is well known to describe the precession of a spinning particle in external electromagnetic 
fields E and B, as has already been demonstrated by Thomas ||27|| . In |E3| it is rederived 
in a manifestly covariant formulation, and therefore also goes under the notion of BMT 
equation. In the present context, ( 4.171) can be viewed as a manifestation of the Ehrenfest 



theorem for spin because here s stands for any expectation value of the spin operator E. 
Notice that since we deliberately omitted a factor of h/2 in the definition (|4.13|) of the spin 
operator, the classical spin vector s is of (constant) unit length. 

So far we have provided geometrical and physical interpretations of d+. Yet, for the 
purpose of the semiclassical representation ( |3.27| ) of the time evolution kernel one needs to 
calculate the SU(2) matrix 

d + (P(t),X(t)) =: (j® , \a(t)\ 2 + \(3(t)\ 2 = 1 , (4-18) 

for each classical trajectory 7^. This can be achieved along the lines presented above, if 

— * — * 

one chooses the initial condition u(£j,y) = ei(£j,y) for a section u solving V + u = 0. In 
the two-spinor representation for u this initial condition reads (wi,o, ^2,0) = (1, 0), and thus 
(O) implies that 



a(t)\ = ( Ul (P(t),X(t)) 

mj \u 2 (P(t),x(t)) 



(4.19) 
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The spin expectation value s(t) corresponding to this choice then follows from (|4.16|) as 



, m x /2Re(a(t)(3(t))\ M 
s(t) = (a(t),P(t)) * R (l = 2 Im («(*) ^W) > with ? (°) = • ( 4 " 2 °) 

w VK*)I 2 - l/W/ W 

Mathematically, this physically motivated construction can be viewed as a map from d + G 
SU(2) to s G S 2 C M 3 . Indeed, it is known as the Hopf map iin '■ SU(2) — > S 2 , which yields 
a U(l) principal fibre bundle over S 2 . 

In order to calculate d + (t) for a given trajectory 7+ y one first has to integrate the evo- 
lution equation (|4.17 ) for the precession of a classical spin along 7+ with initial condition 



s(0) = (0,0, 1) T . According to the relation s(t) = iiH{d+{t)) one thus has determined 
two of the three real degrees of freedom of d + (t) by classical means. The third degree 
of freedom can only be reconstructed, if one is able to lift the curve s(t) in S 2 to SU(2). 
To this end one requires a connection on the U(l) bundle over S 2 that is provided by the 
Hopf map. The two degrees of freedom of the classical spin s can be related to convenient 
coordinates on S 2 once one notices that the normalisation forces |a| to range in the interval 
[0, 1]. Therefore, a variable 9 G [0, it] can be introduced such that 

|a| = cosf and |/3| = sin§. (4.21) 

Together with the phases of a = |a|e ir? and (3 = \/3\ e lA one thus has three degrees of 
freedom at hand to describe d + . A representation of the classical spin s in terms of the 
variables (9,rj, A) then follows from ( |4.20|) , 

'sin 9 cos< 

sin 9 sin (f> ] , with <p := A — r) , (4.22) 
cos 9 

so that (9,(f>) can be identified as the usual spherical coordinates on S 2 . In principle, one 
could now choose A + ^ as the additional spin degree of freedom discussed above. However, 
in order to identify the connection on the Hopf fibration it turns out that a choice of 
either 77 or A separately is more approppriate. Since the north pole of S 2 corresponds to 
(ax, fix) = (1,0) and the south pole to (cts,Ps) — (0, 1), the phase 77 is ill-defined at the 
south pole, whereas A is ill-defined at the north pole. Following the procedure well known 
from the analysis of magnetic monopoles by Wu and Yang , we now choose r\ as a fibre 



coordinate for base points on the northern hemisphere Un C S and, correspondingly, A for 
base points on the southern hemisphere Us C S 2 . Since the connection we want to identify 
is fixed by the requirement that the lifted curve on SU(2) be the solution d + (t) of ( |3.22j ), 
an explicit expression for the connection coefficients can be derived from the equation 



22 



that follows from the identification (|4.19|) . Multiplying this with (ac(t), (3(t)) from the left 
now yields 



Ht)\ 2 - V (t) 



(i 



H*)| 2 ) ^A(t) 



R(t)s(t) 



(4.24) 



Exploiting ( |4.21| ) and (|4.22|) then results in equations for 7/ on Un and A on Us, respectively, 



d 



1 



1 



COS0(t)) -4>(t) 



dt r ] (t) = --R m t)--(l 
^A(t) = -i J*(f)*(f ) + ~ (1 + cos 0(t)) ±6(t) . 



(4.25) 
(4.26) 



dt -w 2 2 r ^" dt 

According to the initial condition in ([4.20 ) the motion of the classical spin s starts at the 
north pole so that for sufficiently small times s(t) G U^. Thus the phase rj should be used, 
whose initial condition 7/(0) = follows from ( |3.21| ). This allows an immediate integration 
of Q, 



T)(t) 



~\ f R{t'W) dt' cos W) ^(0 dt' . (4.27) 



The first term on the r.h.s. of (|4.27| ) is a dynamical phase associated with the classical 
energy of a magnetic moment in given electromagnetic fields, whereas the second term[| is 
a geometric phase. The latter can be further characterised once one takes into account that 
7/ is ill-defined at the south pole so that the phase A should be used instead of rj as soon 
as s enters Us, say, at a time to- One then has to integrate (|4.26| ) with initial condition 
A(to) = <f>{to) + v(to)- To this end the one-form — 1(1 — cos6')d0, which constitutes the 
geometric part of ( [4.27] ), has to be replaced by |(1 + cos9)d(p, see ( |4.26| ). Since these two 
expressions are the gauge potentials of a magnetic monopole, see , we can now formally 
identify the geometric part of ( |4.27| ) to be caused by a magnetic monopole of strength —1/2 
situated at the origin of the sphere. We again emphasize that in the above consideration 
no adiabaticity assumption was made. Therefore, although the result for the geometric 
phase is strikingly similar to the Berry phase of a precessing quantum mechanical spin 
|2T| , the geometric part of r] is rather of the more general type discussed by Aharonov and 



Anandan [|32 



Having integrated the equations for the classical spin and for the additional phase, we 
are now able to present the SU(2) matrix d + in the form 



d + (P(t),X(t)) 



cos 
sin 



(!) 

(f)e 



Kv+<t>) 



sin 

cos 



(I) 
(I) 



-hi 



(4.28) 



where (6,4>) are spherical coordinates for s(t), and i] is given by ( 4.271 ), if s(t) G U^, and 
by f]{t) = A(t) — 4>{t), if s(t) G Us- With this explicit formula for d + all terms entering 
the semiclassical propagator (|3.27|) to leading order in H are completely defined in terms 
of classical quantities. 

1 We remark that in contrast to equation (21) of our recent Letter [Q the sign of this term should be 
as in (Egg). 
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5 Semiclassical trace formula 



After having obtained the fairly explicit expression (|3.27 ) for the time evolution kernel 
together with an interpretation in terms of classical quantities, our ultimate goal now is to 
set up a semiclassical trace formula for the Dirac equation. In the case of the Schrodinger 
equation, Gutzwiller's original approach was to express the quantum mechanical density of 
states in terms of a sum over classical periodic orbits (or, more generally, over connected 
manifolds of periodic points of the classical flow). To this end he Fourier-transformed the 
semiclassical expression for the retarded time evolution kernel, which he derived from its 
path integral representation, in order to obtain a semiclassical approximation for the Green 
function. Subsequently he performed the trace integral with the method of stationary 
phase. The result then immediately yielded the semiclassical spectral density, see [10| for 
details. Since that way one has to deal with several singular objects, we here prefer to derive 
a regularised trace formula that only takes finite quantities into account. The key relation 
for this procedure is equation (|2.27|) , in which we below use a semiclassical representation 
for the truncated time evolution kernel K x (x,y,t). We prefer this procedure, since first of 
all the energy localisation introduced through the truncation x ensures that only the point 
spectrum of H D contributes. Secondly, the test function g cuts off all periods of classical 
periodic orbits that exceed some maximal value T max , since we request g to be compactly 
supported. This cut-off prohibits any possible clash of the two non-commuting asymptotics 
h ^ and t — > oo. Later g can be chosen to have an arbitrarily large, though compact, 
support. Possibly, this support condition can be weakened by demanding a sufficiently 
strong decrease of g{t) as \t\ — > oo, see [[45], pl| . 



The semiclassical analysis presented at the end of section y already revealed that, 
as h — > 0, all contributions to (|2.27|) which exceed 0(h°°) have to derive from classical 
periodic orbits 7^. These are associated with the stationary points (£ 7 ± , x 7 ± , T 7 ± ) of the 
phase (j)^ + Et appearing in the integral (|2.28|) . We also pointed out that the manifolds 
Mq of stationary points (£0, ^o, 0) are related to the hypersurfaces of constant energy, 
which are composed of periodic points with trivial periods T = 0. We now assume the 
cleanness condition explained in section ^ and recall that in particular this implies that 
all further stationary points are such that the periods T 7 ± do not accumulate at t — 0. 
It is then possible to choose a smooth function h G Cq°(R) whose (connected) support 
contains the period To = 0, but no further period T 7 ± > 0. We furthermore require h to 
fulfill h(t) = 1 on some neighbourhood of T = 0. Upon introducing the partition of unity 
1 = h(t) + [1 — h(t)] under the integral ( |2.27| ), the stationary points (£ , ^o, 0) are separated 
from the further stationary points (£ 7 ± , x 7 ± , T 7 ± ) associated with the non-trivial periodic 
orbits •jp. These two classes of stationary points will contribute to (|2.27 ) in essentially 
different ways so that we will deal with them separately. Let our first concern hence be 
the calculation of the leading order contribution of the stationary points (£ ,x ,0). We 
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therefore consider the integral 



h{t) g(t) tr \ X (H D (V £ <P + ,x)) a+ e^ ++Bf ) 



+X [H D (V 2 



(5.1) 



dtd 3 £d 3 x . 



Since the cut-off function h ensures that no stationary points with t 7^ contribute, we 
expand the phases about t = 0. Using the Hamilton- Jacobi equation for S ± (x, £,£), one 
finds 



(x, x, t; £) + Et = S ± (x, f , t) + Et- x£ 
= t 



dS ± - ' 
E + —(x,£,0) 



Of 



+2 ff. q± 

+ 2^r(£,£,o) + o(* 3 ) 



(5.2) 



t 



f 2 



E - , x)\ + -VpH^Z, x) V 3 H±& x) + 0(t 3 ) 



At this stage we introduce in ( J5.1[ ) polar coordinates for the variable £, i. e., £ = Au> with 
A := |£| and solid angle w e S 2 C 1R 3 . This implies d 3 £ = \ 2 d\duj. We then employ the 
method of stationary phase to the integration over the variables (t, A) G M x R + . Stationary 
points are therefore determined by a vanishing of the derivatives of ( |5.2|) with respect to 
t and A, respectively. Evaluated at t = 0, the expression ( |5.2j ) shows that this yields the 
condition if ± (Ag cu, x) = to be fulfilled by the stationary points (tj = 0, Aq ). Therefore, 
our first conclusion is that the integral over the remaining variables (x,u) 6 K 3 x S 2 is 
in fact restricted to the hypersurfaces fl E and fl E , respectively. Carrying out the method 
of stationary phase further finally yields an asymptotic expansion of the integral ( |5.1|) as 
h ^ whose leading term can be determined explicitly in a straight forward manner, 



2tt 



where vol(fi^) denotes the volumes, 



0(0) 2vol (g|) + 2vol (Q E ) 



0(h)} 



(5.3) 



vol(f2|) = I [ 5(H ± (p,x) — E) d 3 pd 3 x 

./R 3 ./R 3 



(5.4) 



of the hypersurfaces 

For the computation of the contribution to (|2.27|) caused by the second class of station- 
ary points, associated with the non-trivial classical periodic orbits, we proceed differently. 
As a starting point we consider the r.h.s. of ( 2.27|) , in which we introduce the semiclassi- 
cal representation ( |3.27| ) of the truncated time evolution kernel. That is, we are going to 
evaluate 



1 

2^ 



+00 



[1 - h{t)} g(t)e* Et UK x (x,x,t) dtd 6 x 



(5.5) 
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with the method of stationary phase, after having inserted the r.h.s. of ( |3.27| ) for K x (x, x, t). 
Notice that the factor 1 — h(t) cuts off a neighbourhood of t = and that the compactly 
supported function g further restricts the range of integration over t to a bounded set. 
The semiclassical asymptotics ( |3.27|) , which was derived with t held fixed, can therefore 



be used in ( p.5[ ) without running into conflict with a need to perform t — > oo. Apart from 
the factors caused by the spin degrees of freedom in ( |3.27| ), and from the fact that the 
two relativistic classical Hamiltonians H ± (p, x) determine the equations of motion for the 
translational degrees of freedom, the following calculation is closely parallel to the case of 
the Schrodinger equation flO| . 



As implied by the discussion at the end of section 0, the stationary points relevant for 
the semiclassical evaluation of (|5.5|) , which we now express as 



27r(27rift) 



J-oo + 

-)tr (w t d„Wl) £T_ e i[R ^ +m ~'^ x 1 {1 + 0(h)} dt d 3 x , 

xx \ / 'Jxx I 



(5.6) 



derive from the non-trivial periodic orbits 7^. Indeed, the stationary points of the phases 
appearing in (|5.6|) are determined by 



VgR± (x, y, t) + VgR± (x, y,t) =0 and (x, x,t) = E . (5.7) 



y=x 



dt 



The first relation requires that initial and final momenta of those closed trajectories 7^. 
that contribute to the method of stationary phase must coincide, so that these are indeed 
periodic. The second condition then picks out those periodic orbits that are contained in 
the hypersurfaces At stationary points the phases then read 

R % (V . % > T ,t ) + E T ,t =<fj&2=-- < ± (E) ■ (5-8) 

In order to apply the method of stationary phase to fl5.6| ) we assume that all periodic orbits 
are isolated and non-degenerate. Strictly speaking, one now has to introduce a partition of 
unity that separates the contributions of the isolated orbits. Since their periods are known 
not to accumulate at some finite value T, this can be done in the same manner as for the 
trivial period Tq = above. 

The only difference to the case of the Schrodinger equation is caused by the presence 
of the spin degrees of freedom. This is represented through the matrix trace under the 
integral in ( |5.6| ). At stationary points, however, t = T 7 ± and the periodicity implies that 
Vt , — Vq. Thus a cyclic permutation under the matrix trace together with the former 
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result qPSP yields 

tr(Vr + d +at Vj) = tr(yj V d + ^) = tr (d + , 7 +) = 2 cos(^ >+ ) cos V . (5.9) 



Finishing the calculation as in the well known case of the Schrodinger equation || |TU 
finally yields a contribution of 

X(E) A yt e^ iE) {1 + 0(h)} (5.10) 

for every isolated non-degenerate periodic orbit 7^. Here the amplitude A^± contains only 
classical information about the periodic orbit, including the contribution of a classical spin 
precessing along the orbit. Explicitly, the amplitude reads 



2T\ cos(^# +) cos 77 ± 
A 7 ± = ^^e >. (5.11) 



det(M >± - 1 



1/2 



In this expression T* denotes the primitive period of 7^, i. e., the period of the associ- 
ated primitive periodic orbit. Moreover, M^± is the linearised Poincare map (monodromy 
matrix) along the orbit and /i 7 ± is its Maslov index. 

From the Gutzwiller trace formula it is well known that if 7^ is not primitive, all 
quantities that enter ( |5.10| ), apart from the spin contribution ( |5.9|) , can be readily expressed 



in terms of the respective quantities of the associated primitive periodic orbit. In order to 
extend this to ( |5.9| ) we recall that if 7^ is a fc-fold repetition of a primitive orbit, where 
fc6Z \ {0}, the fact that d +1 + is a holonomy implies that d +1 + = {d^+) k , where ^* 7 + 



denotes the corresponding primitive holonomy. According to [46] the trace Q5.9Q associated 
with 7^ can hence be expressed as 



where T^{x) is the fc-th Chebyshev polynomial of the second kind in x. Thus, if desired, 
the following trace formulae ( |5.13| ) and ( |5.15| ) can also be expressed in terms of sums over 
primitive periodic orbits and their repetitions. 

In case the classical dynamics generated by the Hamiltonians H ± have only isolated 
non-degenerate periodic orbits, the relations (|5.3| ) and ( p.lOj) now enable us to state the 
following semiclassical trace formula explicitly, 

E n - E\ _ |(0) 2 vol (g|) + 2 vol {yt~ E 
X[ ' 2tt (2tt/i) 2 



2tt > 

7 P 
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The conditions imposed on the test function q and its Fourier transform ensure that all 
expressions entering this trace formula are finite. In particular, due to the compact support 
of q the sum over periodic orbits only includes orbits up to a finite period and since the 
Hamiltonian flow generated by H (p,x) was supposed to be clean the sum is of finite 
length. 

Often a semiclassical trace formula is presented for the spectral density of the quantum 
Hamiltonian, see e. g. H |i~0fl . In the present case one can readily obtain such a trace 
formula for the truncated spectral density 



d x (E) :=J2x(En)S(E-E n ) (5.14) 



from ( |5.13| ). This reads 



1 ^ T | cos (|^)cosr/ 7 ± ks±{E) ^ ± ^ (5-15) 



± |det(M >± -l 4 ) 



Obviously, the sum over classical periodic orbits in ( |5.15|) does not converge. This trace for- 
mula rather has to be viewed as a distributional identity whose actual meaning is provided 

by (E13D . 



6 Nonrelativistic limit 

As compared to the Schrodinger equation, the Dirac equation takes two generalisations 
into account. It first takes care of the spin degrees of freedom and, secondly, describes 
relativistic dynamics. In many situations of physical interest it, however, suffices only to 
include spin and to leave the description of the translational motion on a nonrelativistic 
level. As it is well known, this can be achieved by considering the Pauli equation, possibly 
with several additional terms such as one describing spin-orbit coupling. In this section we 
therefore want to study particles of charge e and mass m with spin 1/2 in a nonrelativistic 
context. To this end we compare the semiclassical asymptotics for the Pauli equation with 
the leading order as c — >• oo of the time evolution kernel ( p.27| ) for the Dirac equation. 
Both approaches, which will turn out to produce identical results, then allow to set up a 
semiclassical trace formula. 

We recall that when one divides a Dirac four-spinor \I/ into two-spinors according 

to 

•<*•«>-= fell!) • <"> 
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one obtains from ( |2.1Q two coupled matrix differential equations. It is well known []36 
that to leading order as c — > oo these equations decouple, and that for ip := ip 1 the Pauli 
equation 



eh 
2mc 



oB(x) 



ip(x,t) =: H P ip(x,t) (6.2) 



emerges. See also |38l for a careful treatment of the limit c — > oo. In Weyl quantisation 
the Pauli Hamiltonian Hp can be realised as 



H, 



Hi- I -V^, X 



with the 2x2 symbol matrix 

H P (p, x) = H (p, x) + hHi(p, x) 



(6.3) 



(6.4) 



As opposed to the Dirac equation, compare (|2.5|) , this Weyl symbol is composed of a 
principal symbol 



H (p,x) :-- 



^(p-- c m) 2 + e<p(x) 



1 



2 j 



which is a multiple of the identity matrix, and an additional subprincipal symbol 



Hi(p,x) 



2mc 



aB(x) 



(6.5) 



(6.6) 



that reflects a coupling of the spin degrees of freedom to the external magnetic field. 

We are now going to study the semiclassical limit h ^ along the lines developed 
above for the Dirac equation. Since therefore many details of the following calculations are 
similar to the ones shown in the previous sections, at several places the presentation will 
be kept rather brief. The matrix- valued Schwartz kernel Kp(x,y,t) of the time evolution 
operator Up{t) := e~h Hpt obeys 



K P (x, y, t) V>o(y) d 3 y , ip(x, 0) = ip (x) 



ip(x, t) 



so that Kp(x,y,t) has to solve the Pauli equation (|6.2|) with initial condition 

4 hm K P (x, y, t) = 1 2 6(x - y) . 



(6.7) 



(6.{ 



Since the principal symbol Ho is a scalar multiple of the identity matrix and thus has one 
eigenvalue, only one Hamilton- Jacobi equation will be relevant to the semiclassical time 
evolution. We therefore choose the semiclassical ansatz 



K P (2, y, t) = -i— / f; (?) a k (x, y, t; f) e*^)"^ 
(2tt^) 3 J r3 V 1 / 



d 3 ^ + 0{h c 



(6.9) 
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compare ( |2.8| ) and ([2.14J) , which we introduce into ( |6.2| ). Comparing like orders in h yields 
to lowest order a Hamilton- Jacobi equation for S, whereas to orders h k , k = 1,2,3, ... , 
transport equations for the 2x2 matrices a>k-i follow. The initial conditions for S — y£ 
and dk are analogous to ( |2.10|) and CP .11 ). By obvious reasons, the subprincipal symbol 
( |6.6D cannot appear to zeroth order in h so that the Hamilton- Jacobi equation emerging 
in leading semiclassical order is only determined by the principal symbol H (p,x). As a 
consequence, the classical dynamics of the translational degrees of freedom are those of 
a nonrelativistic point particle that does not experience a force coming from a coupling 
of spin to the external magnetic field. The latter is, however, contained in the transport 
equation for a that appears in next-to-leading order. This reads 



d 1 - 

V P H (V?S, x)V s + — + - V* V f H (V 3 S, x 



+ iMp(x) 



a Q = 



(6.10) 



so that the spin degrees of freedom enter through the traceless hermitian 2x2 matrix 

e 



Mp{x) :-- 



2mc 



aB(x) = H^V^x) . 



(6.11) 



In analogy to the considerations following ( |3.14| ), the ansatz 



det 



(■ 



d 2 S 



\dx k d£i 



(£,Z j ,t))dp(P(t),X(t)) 



(6.12) 



for the lowest order amplitude, evaluated along a classical trajectory ^ xy associated with 
the stationary point £j of the phase, proves useful. It leads to the nonrelativistic spin 
transport equation 



- + xM P {X{t)) 



d P {P{t),X{t)) = , d P (P(0),X(0)) = 1 



2 • 



(6.13) 



along j xy that determines the SU(2) matrix dp nxy describing the leading contribution of 
spin to the semiclassical time evolution kernel. The transport equations Q6.10 ) and ( |6. 13 ) 
have already been obtained by Choquard J47|, who used a semiclassical ansatz similar to 
( |B~9|) , however, without the integration over £. 

We remark that the decomposition of the SU(2) matrix d + into a classical spin s and 
an additional phase rj presented in section [| can be repeated in the present, nonrelativistic, 
context. This procedure leads to the dynamical equation 



dt 



sx — B 

mc 



(6.14) 



for s transported along j xy in the external magnetic field. The r.h.s. of (|6.14j) can readily 
be identified as the leading order of the r.h.s. of (|4.17|) as c — > oo. An additional phase, 



vp(t) 



2mc 



B(t')s{t') dt' 



(1 



C0S e(t'))^(t')dt' , (6.i5) 
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representing the third, nonclassical, degree of freedom of dp, appears in the same manner 
as for the Dirac equation. This phase is again composed of a dynamical and a geometric 
part. 

A calculation similar to that presented in section |^ finally leads to the following semi- 
classical time evolution kernel for the Pauli equation (for t ^ 0), 



K P (x,y,t) 



1 



(27ri£) 3 / 2 



5> 



Ixy ^fxy 



n v ^v{i + o{h)} 



(6.16) 



Ixy 



On the r.h.s. the sum extends over the solutions of Hamilton's equations of motion with 
the classical Hamiltonian H (p,x), which connect y and x in time t. The corresponding 
Morse indices and Hamilton's principal functions are denoted by v lxy 
respectively; compare (|3.18j ). The factor _D 7 is defined by 



and R lxy (x,y,t), 



D 



Ixy 



det 



d 2 R 



-y X y 



dx k dyi 



x, y, t) 



1/2 



(6.17) 



Since these quantities are already determined by the nonrelativistic classical dynamics 
generated by H , the only difference between (|6.16|) and the respective semiclassical kernel 
for the Schrodinger equation is an appearence of the factors dp >lxy G SU(2) that represent 
the leading influence of spin. 

So far we have examined the propagator that emerges from first taking the nonrela- 
tivistic limit of the Dirac equation and then constructing a semiclassical time evolution 
kernel. We will now compare this to the result that one obtains by first constructing the 
relativistic semiclassical kernel ( pj.27| ) and then taking the nonrelativistic limit. We are 
thus now interested in the leading order behaviour of ( |3.27| ) as c — > oo. To this end we 
first remark that 



V(p, x) 



(o)+°0 , W(p>-)=(°)+0(i) 



oo 



(6.18) 



i. e., in leading order the positive and negative kinetic energy parts in (|3.27|) decouple 
completely. We therefore have to compare ( |6.16| ) with the upper left 2x2 block of (|3.27| ) . 
This is consistent with (|6.1| ) because this block describes the time evolution of the two- 
spinor ip = ipi. We furthermore recall the well known fact that the classical relativistic 
dynamics described by H + (p,x) turns into the nonrelativistic dynamics with the classical 
Hamiltonian H Q (p,x) as c — » oo. This in particular implies that 7+ R + + , z/ + + and D + + 

Ixy Ixy 



may in leading order be approximated by 7^, R 1 , v^ xy and D lx . We are therefore 
only left with comparing the factors d + and dp which describe the influence of spin in 
a relativistic and in a nonrelativistic context, respectively. As already mentioned below 
equation ( |6.14| ) 



Mi 



M P + oa) 



00 



(6.19) 
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leads to a nonrelativistic approximation for both the classical spin dynamics and the non- 
classical phase. Thus dp provides the leading order asymptotical term for d + as c — > oo. 
Collecting everything one observes that in the nonrelativistic limit the upper left block 
of the semiclassical time evolution kernel ( 3.27Q for the Dirac equation turns into the re- 



spective result ( |6.16| ) for the Pauli equation. In this sense the limits h — > and c —>■ oo 
commute, at least concerning leading orders. 

In order to set up a semiclassical trace formula for the Pauli Hamiltonian we now 
assume that Hp has a pure point spectrum. Otherwise we would have to employ an 
energy localisation to a gap in the essential spectrum as described in section |2| for the 
Dirac Hamiltonian. We then consider a test function g £ <S(M) with compactly supported 
Fourier transform g. Applying the procedure of section |5| to Hp finally yields 

V- ( E n -E \ e(0) 2vol(n E ) M 

^ 6 {—n-) = l27nF { + ( )} 

(6.20) 

+ £^J^ , e *^{l + 0(fi)}, 

z — ' Z7T 

where the sum extends over the periodic orbits of the classical dynamics generated by the 
principal symbol Hq. In fact, all classical quantities entering ( |6.20| ) refer to this Hamil- 
tonian. As explained in the relativistic case, for such a trace formula to be valid the 
flow generated by Hq must be clean. The amplitude associated with each isolated, non- 
degenerate periodic orbit then reads 

2Tf cos(i6> ) C0S77p~ 
= _> "lUll J£l1l e -if^ . (6.21) 

|det(M 7p -l 4 )| 1/2 

A different semiclassical approach to the Pauli equation has previously been used to 



investigate spin-orbit coupling [^, 34]. The authors of these papers principally base their 
method on the technique developed by Littlejohn and Flynn || [7]] to treat matrix-valued 
wave operators with principal symbols that have no (globally) degenerate eigenvalues. 
Below we will derive the time evolution according to the prescription of the semiclassical 



limit employed in |34|, however, using the techniques outlined in the previous sections. 
For simplicity, and for ease of comparison with our previous semiclassical study of the Pauli 
equation, we will not consider spin-orbit coupling but only a coupling of spin to the external 
magnetic field. Therefore, the relevant Hamiltonian is the one defined in fl6.2|) . However, 
a generalisation of the following discussion to arbitrary couplings of the form aC(p,x) is 
straight forward. Following now the philosophy of |33|, [34]], we introduce Bohr's magneton 
a := anc i consider it as constant in the semiclassical limit. Thus the Hamiltonian 

~ Zmc 

1 (h _ e - \ 2 
H' P := — -V - -A(x) + e <p(x) - u BBix) (6.22) 
2m \ l c / 

arises as a Weyl operator associated with the symbol 

H' P (p,x) := — (p- -A(x)) 2 + e<f(x) - /iaB(x) , (6.23) 
2m \ c J 
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so that no subprincipal symbol occurs. As opposed to the situation analysed at the be- 
ginning of this section, one could view the present procedure as taking the simultaneous 
limits h ^ and \B\ — ► oo in such a way that H\B\ — const. Another way to look at this 
is to keep \B\ fixed, but to perform the limit of 'large spin'. 

As long as B ^ the 2x2 symbol matrix ( p.23| ) has two non-degenerate eigenvalues 



H'%{p,x) 



2m V c 



+ e<p(x) =F /i \B(x) 



(6.24) 



which generate two classes of classical dynamics. Since this is similar to the situation 
occurring for the Dirac equation, the semiclassical ansatz for the time evolution kernel 
K' P (x, y, t) therefore should be chosen as in ( [2.8| ), but where now the amplitudes are 2x2 
matrices. Applying the procedure described in section |3| then leads to two Hamilton- Jacobi 
equations, with the two Hamiltonians ( |6.24| ). 

The non-scalar contribution to the symbol matrix (|6.23| ) is given by the Hamiltonian 
— n aB(x) describing a quantum mechanical precessing spin. The ortho normal eigenvectors 
v±(x) G C 2 of H'p(p,x) corresponding to the eigenvalues H'^,(p,x) are hence well known 
from the standard example of the Berry phase [ 29fl . In analogy to (|3.5| ) one can now 
introduce the ansatz 



(6.25) 



with some suitable vector v G C 2 . Upon multiplication of the transport equation for 
with v+(x) from the left one obtains the equation 



V P ^(V^ + , x)V 3 +- + [V f H' P (V £ S- 







Tl+I 



,x 



iM'+(V^ + ,x) 



5+ 



, (6.26) 



where 



M'p(p } x) := -— (p- -A(x)) vUxWgvJx) . (6.27) 
m V c / 

On the r.h.s., the factor ulVv+ is the well known expression for the adiabatic connection 



leading to the U(l)-Berry phase of a precessing spin ||29|1 . Again, the transport equations 
for will be solved along the classical trajectories 7' following from the Hamiltonians 
( |6.24j ). As suggested by (|3.15|) and ( |3.20|) , one separates 



vt(j£(t),$,,t) = Jdet (J^^Xit),^ d' + (P(t),X(t)) vl(y) , (6.28) 



so that along a trajectory 7'^ the phase d'+(t) G U(l) has to solve the equation 



- + vi(X(t))V,v + (X(t))X(t) 



d'+(t) = , d'+(0) = 1 . 



(6.29) 
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As a result, one obtains the well known U(l)-Berry phase of a quantum mechanical spin 

/+ 

xy 



that has been transported adiabatically along the trajectory 7'™. Finally, the semiclassical 



time evolution kernel reads 



K ' p[i ' l) = (27rife)3/2 E D y± e*^" 1 *"^ {1 + 0(h)} , (6.30) 
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where all classical quantities refer to the dynamics generated by the Hamiltonians ( |6.24| ) 



and are defined in analogy to the previous cases. In particular, Hamilton's principal func- 
tions R'y± (x, y, t) are related to the analogous quantities Ry±^(x, y, t) that are defined by 

the Hamiltonian fl6.5|) , but here evaluated along the trajectories 7' through 



R*±{x,y,t) = R 1% (x,y,t)±fi f \B{X{t'))\ dt' . (6.31) 

J 

Due to the factors v±(x)v±(y) in (|6.30|) the 'spin up' and 'spin down' components of 
an initial spinor ipoiv) are propagated independently along the corresponding trajectories 
following from the Hamiltonians if'p. Here 'spin up' and 'spin down' are defined with 
respect to the instantaneous direction of the magnetic field. This procedure breaks down 
at mode conversion points, i. e., at points where the magnetic field vanishes. There the 
symbol matrix ( |6.23j ) has one twofold degenerate eigenvalue and the level surfaces of the 
two Hamiltonians H' p cross. Moreover, since V^-fTp develops a singularity, the classical 
trajectories are not smooth when crossing a mode conversion point. In certain situations 
this defect can be cured by letting the trajectories cross the two level surfaces, see [|34| , 
but in general an application of the present semiclassical procedure requires a more refined 
treatment of mode conversion points, see for a detailed discussion. 

At this place one could easily establish the corresponding semiclassical trace formula, if 
one followed the programme outlined in section |2] once again. However, we refrain from do- 
ing this here and rather comment on the relation between the two semiclassical approaches 
to the Pauli equation discussed in this chapter, which lead to the two distinct expressions 
( |6.16| ) and ( |6.30| ) for the time evolution kernel. In the first scenario we systematically per- 



formed an expansion in h and determined the leading order terms for the time evolution 
kernel and for the classical side of the trace formula. We observed that to lowest order 
the translational degrees of freedom decouple from the spin degrees of freedom in that the 
translational motion experiences no back reaction from the coupling of spin to the external 
magnetic field. In a certain sense this decoupling can be seen as an adiabatic one where 
the translational motion is considered as slow, although this condition is not needed for 
the formulae to be valid. This topic was also discussed by Balian and Bloch p|J in the con- 



text of semiclassical approximations for the Green function. Spin enters in next-to-leading 
order and, among other quantities, determines the amplitudes in the relevant semiclassical 
expressions. The leading order of the spin dynamics is given by that of a classical spin pre- 
cessing along the particle trajectories. We have repeatedly emphasized that at this stage 
no adiabatic limit is considered. In addition, a geometric phase of the type discussed by 
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Aharonov and Anandan |32] enters the amplitudes. In the second scenario we considered 
the double limit h — *■ 0, \B\ — > oo, ft|.B| = const. It turned out that in this context the 
relevant classical translational motion follows from two Hamiltonians, taking the effect of 
a coupling of a 'spin up' and a 'spin down', respectively, to the external magnetic field into 
account. Since in this context, via the expressions v±(x)v±(y) in (|6.30|), the spin direction 



is defined with respect to the instantaneous direction of the magnetic field, the spin degrees 
of freedom are transported adiabatically along the particle trajectories. There is no further 
dynamical equation for a classical spin and, as a consequence of the limit \B\ — > oo, the 
Berry phase that enters can be viewed as emerging from an adiabatic approximation of the 
geometric term found in the first scenario. This finding is in agreement with the remark 
on the relation of these phases that can be found in ||32|| . In order to be more specific 
concerning the issue of adiabaticity, we introduce spherical coordinates ($b,¥b) for B/\B\ 
which allows to obtain an explicit expression for the Berry phase emerging from (|6.29| ). 



Moreover, we reintroduce for the magneton \i. Using (|6.31| ) one thus observes that 

d' ±e * R '^« = exp(±i— I |B|dt'=F- / (1 - costf B ) u> B dt'l e*M, , (6.32) 
I 2mc J 2 J J 

which can be introduced in ( |6.30| ). The first factor on the r.h.s. is readily identified to arise 
as an adiabatic approximation to e ir,p , if the phase (|6.15|) is evaluated along the trajectory 
"f'xy instead of ^ xy . Thus the spin contributions to ( |6.30|) are adiabatic approximations to 



the respective contributions to (|6.16|) . Conversely, one concludes that in the semiclassical 



time evolution according to ( |6.16| ) the spin transport is performed non-adiabatically so 
that mode conversion poses no difficulty In conclusion, one can consider the two alter- 
native ways of performing the semiclassical limit as being considered with a 'weak' and a 
'strong' coupling, respectively, of spin to the translational degrees of freedom. Comparing 
the two results fl6.16| ) and ( |6.3(J| ) for the time evolution kernel one notices that the two 
approaches yield different results when extended to intermediate couplings. In the first 
case the translational motion is not influenced by the spin, whose dynamics in turn is not 
treated adiabatically. In the second case, however, spin has an effect on the translational 
motion, but the spin dynamics enters in an adiabatic approximation. This observation 
illustrates the fact that the ^-expansions employed are not uniform in the field strength. 
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